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Truncated sum rules have been used to calculate the fundamental limits of the nonlinear suscep- 
tibilities; and, the results have been consistent with all measured molecules. However, given that 
finite-state models result in inconsistencies in the sum rules, it is not clear why the method works. In 
this paper, the assumptions inherent in the truncation process are discussed and arguments based on 
physical grounds are presented in support of using truncated sum rules in calculating fundamental 
limits. The clipped harmonic oscillator is used as an illustration of how the validity of truncation 
can be tested; and, several limiting cases are discussed as examples of the nuances inherent in the 
method. 
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I. INTRODUCTION 

Truncated sum rules were first used at the turn of 
the century to calculate the fundamental limits of the 
ofi- resonant hyperpolarizabiUty /3,P, 0] and the second 
hyperpolarizalDility, 7.(11 0| The role of off-diagonal el- 
ements, such as measured with Hyper Rayleigh Scat- 
tering, were also considered. These calculations were 
later applied to calculating the resonant two-photon ab- 
sorption cross-section 6| as well as the maximum possible 
two-photon absorption cross-section when a molecule is 
doubly-resonant with the excitation source. 7] The the- 
ory is supported by the fact that no molecule has ever 
been found to break the fundamental limit. Q 

While all molecules ever measured were found to fall 
below the limit, it was pointed out that there is a gap be- 
tween the fundamental limit and the molecules with the 
largest measured nonlinear-optical susceptibility. |^ ^ |^ 
This apparent limit, a factor of 10'^/^ below the funda- 
mental limit, had two possible implications. Either, the 
theory could be flawed so that it overestimated the fun- 
damental limit by more than a factor of thirty; or there 
was another factor that acted to suppress the nonlinear- 
optical response. To test the hypothesis that the theory 
has overestimated the nonlinear response, all one needs to 
do is show that a system exists whose hyperpolarizabiUty 
exceeds the apparent limit. The clipped harmonic oscilla- 
tor, whose hyperpolarizability can be calculated analyt- 
ically without approximation, has a hyperpolarizability 
that is an order of magnitude larger than the apparent 
limit, yet about a factor of two below the fundamental 
limit. [IG] Given that very talented organic chemists have 
been working almost three decades to improve the hy- 
perpolarizability of molecules through structure-property 
studies, and no molecule has ever been found to breach 
the apparent limit, it is probably prudent to conclude 
that the best organic molecules fall below the apparent 
limit. 

While nanoengineering methods have been used to in- 
crease the second hyperpolarizability per molecule to 



within a factor of two of the fundamental limit, a 
more careful analysis of this case shows that the interac- 
tions between molecules makes them respond collectively. 
As such, the collection of molecules is acting like a single 
supermolecule where its response is still well below the 
apparent limit. So, the state of affairs with regards to 
molecular hyperpolarizabilities is that all molecules ever 
measured and all analytical calculations of P and 7 ever 
performed fall below the fundamental limit, though some 
of the calculations come close. This set of evidence sup- 
ports the theory of fundamental limits. So, while the 
process of truncating the sum rules and applying them 
to the analysis of the Sum Over States (SOS) expression 
for (3 could lead to pathologies, the results of this process 
seems consistent with observation. 



II. THEORY OF THE FUNDAMENTAL LIMIT 



This section develops the theory 
limits of the hyperpolarizability (the 
izability follows along the same lines 
more messy, so will not be presented 
presented follows the literature, P, 01 
will be focused on the assumptions 
culations and all of the details will 
commentary on their meaning. 

We being from the sum rules. 



of the fundamental 
second hyperpolar- 
but is algebraically 
here). The method 
but, close attention 
underlying the cal- 
be presented with 
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which relate the matrix elements of the position operator, 
x„,„, to the energy differences, Enm = En — Em, between 
the states n, and m. Since the sum rules derive directly 
from the Schrodinger Equation, they are of a fundamen- 
tal nature and apply to all systems in which one can write 
a potential energy function. Note that Equation ^ yields 
a distinct equation for each distinct pair of values of m, 
and p. As such, we will refer to the sum rule equation 
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with indices m and p as Equation Xmp- When the sum 
rules are truncated to include only £ terms (i.e. the sum 
over n in Equation granges from to € — 1), we will refer 



to the truncated sum rule as x 



A. Linear Susceptibility 

We begin by treating the linear susceptibility, which 
is the simplest possible case. A quantum perturbation 
calculation of the polarizability of a one-dimensional sys- 
tem, a, yields. 
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(2) 

where w is the frequency of light, ujno the transition fre- 
quency (jJn — ijJq between states n and 0, 7„o the width of 
the transition, and e the electron charge. We consider the 
off- resonance limit, where a; ^ and calculate the real 
part of the polarizability, where i?e[a], where 7„o = 0, 



n—1 



(3) 



where we have used Eno = fu^no- 

Now we use the non-truncated sum rules to simplify 
Equation 01 The sum rule Equation Xqq (Equation ^ 
with m = p = 0) is given by 
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To use Xqo, we rewrite Equation 13 as follows. 



EnO \Xno\ 

E'ip, 



< 



EnO \Xno\ 
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where the inequality follows from the fact that Eno < 
i?io and all the terms in the sum are positive definite. 
Substituting Equation 0] into Equation |3J we get 



Notice that if we write sum rule Equation Xqo as. 



(6) 



l-ior-^iV-El^Kor, (7) 
2mEio Eio 

then clearly, the maximum value of a is obtained by plac- 
ing all of the oscillator strength into the first excited state 
and setting all of the other transition moments to zero, 
that is Xno = for n > 1. Note that when all of the tran- 
sition moment is concentrated into the first excited state, 
xio takes on its maximum value as given by Equation 



The important lesson of this calculation is that in the 
limit of maximal a, the oscillator strength gets concen- 
trated into the transition between the ground and first 
excited state. It is instructive to treat this problem as 
a three-level system where the SOS expression includes 
only the ground and first two excited states; and simi- 
larly, to truncate the sum rules to three levels. 

The three-level model for a is given by 
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where by definition, state 2 is of higher energy than state 

(3) 

1. Using the three- level sum rule given by Equation Xqq 
and including only three states (0, 1, and 2) to eliminate 
xo2, Equation ijSJl becomes: 

2l^iol' Nh^ (Eio\ 1 

Ew \ 2mE2o {e^J \xwf \E2o J J 



Equation (O can be rewritten in the form, 

a3L = aMAX {X^ + E^- E^X^) = aMAxHX, E) 
where 



(9) 
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with 
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and 



E : 



£^10 

£^20 



(11) 



(12) 



(13) 



Note that with the above definitions, — 1 < x < +1 and 
— 1<E<+1. We call Equation [TUI the reduced three- 
level model since the sum rules have reduced the number 
of parameters. 

Using the reduced three-level model, we assume that 
the energy parameter, E, and the transition moment pa- 
rameter, X, are independent variables. Figure ^ shows 
a plot of h{X,E) as a function of X for several values 
of E. Clearly, h{X — l^E) = I for all values of E, so, 
the three-level model yields the same maximum value of 
a = aMAX as the infinite-state result. 

The reduced three-level model yields insights into the 
nature of fundamental limits. First, the fact that X ^ 1 
at the fundamental limit clearly shows that all of the 
oscillator strength must reside in the first excited state. 
This is consistent with the results of the infinite level 
model. Secondly, note that at E = 1, when the two ex- 
cited states are degenerate, the polarizability is at the 
fundamental limit for all possessible values of the tran- 
sition moment to the first excited state. This behavior 
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FIG. 1: A plot of h{X,E). 
can be understood by noting that when Eiq = £'20, the 

(3) 

truncated sum rule Xqq yields 

M' + \x2o\' = :r^^N, (14) 
which when substituted into Equation |H1 yields 

£^10 1^ J mSfo 

As such, the degenerate case acts like a two-level model. 

This simple example shows that the fundamental limit 
of the polarizability occurs when all of the oscillator 
strength is concentrated in the transition to the first ex- 
cited state. As such, it is sufficient to consider only a 
reduced three-level model since this is a good approxi- 
mation for any quantum system that has a polarizability 
near its maximum value. 



B. Nonlinear Polarizability 

In the case of the hyperpolarizability and higher-order 
polarizabilities, the sum-over-states expressions are so 
complicated that it is not possible to calculate the exact 
result as we did for a. In the linear case, we concluded 
that the fundamental limit occurs when all the oscilla- 
tor strength is concentrated into one transition. How- 
ever, the nonlinear susceptibilities depend on two types 
of terms. So, it would not be unreasonable to expect that 
the oscillator strength may need to be shared between at 
least two excited states. 

Assuming that a fundamental limit indeed exists, we 
therefore use the following ansatz:p^ 

• Ansatz: In the limit when the static nonlinear 
polarizability of a bound quantum system in its 
ground state approaches the fundamental limit, the 
only allowable non- vanishing matrix elements of the 
dipole operator are described by a 3x3 matrix. 



An important consequence of this ansatz is that the sum 
rules must be truncated, resulting in a finite number of 

(3) 

equations of the form ATAm, with to < 3 and n < 3. The 
set of truncated sum rules leads to several problems, such 
as self inconsistencies. To deal with inconsistencies and 
decide which set of sum rules are appropriate, we adopt 
the following guiding principles: 

1. The set of sum- rule equations must be self- 
consistent. 

2. The predictions following from the application of 
the sum rules must be consistent with experiment. 

3. When off-resonance, all measurable quantities 
should be free from divergences. 

Our first priority is to decide on the fundamental set 

of sum-rule equations. First, we note that sum rule 

(3) 

equation yields a nonsensical result, and, it is in- 

(3) 

consistent with the sum rule equation Xqq . Based on 

(3) 

Principles 1 and ^, we eliminate . Secondly, we 
note that the remaining sum-rule equations are overly 
restrictive since they predict that all molecules should 
have the same ratio, E. This violates Principle 2 since 
molecules with a continuous spectrum of E vales are ob- 
served. As such, we must drop at least one additional 
equation. To determine which equation to drop, we go 
through the process of dropping one equation, and an- 
alyzing the ramifications of the remaining set. We find 
that the only equation that can be dropped, which leaves 
a set of equations that do not result in divergences in 
the calculated off-resonant nonlinear susceptibilities (i.e. 

obeys Principle 3), is X^^ . 13] If we drop an additional 
equation, then the result one gets is that there is no limit 
on the susceptibility. Since experiment strongly suggests 
that the nonlinear-optical response can not be arbitrar- 
ily large, we conclude that the 4 independent equations 

Aqq\ x'ff^ , X^^ , and X^ form the complete set. (Note 

(3) (3) 

that Equations Xnm and Xmn are complex conjugates of 
each other, and for a real dipole matrix, are identical.) 

The process of applying the truncated sum rules to the 
three-level model of the sum-over-states expression of the 
hyperpolarizability, (3 and second hjnperpolarizability, 7, 
is described in the literature, dimil For^, this yields, 

I ^MAX 1 3 

/? = -e^ ' ^^2^ ' G{X)f{E) = l3oG{X)f{E), (16) 
where 

f(E)^{l-Efl^{E^ + ^E + l^, (17) 

and 



G(A)=^aW|(1-A4). (18) 
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As in the linear case, we assume that E and X are in- 
dependent. f{E) peaks at E = and G{X) peaks at 
X = ~</3. The maximum value of each function is unity, 
so from Equation 1161 we get the fundamental limit 



/3. 



MAX 



: /3o/(0)G( = ^/3 



7V3/2 



7/2 



(19) 



where we have made use of Equation 1121 A similar pro- 
cedure for the second hyperpolarizabihty yields SB 



R5 

^10 



< 7 < 4- 



R5 
^10 



(20) 



Note that the fundamental limit of of positive 7 is four 
times larger than for negative 7. 
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III. ANALYSIS 



FIG. 2; Contour plot of re^p. Inset shows matrix elements. 



In this section, we develop a method for analyzing the 
validity of using truncated sum rules by studying how 
the truncated series compares with the exact results. We 
begin by defining the K-matrix, which derives from the 
sum rules given by Equation ^ 



TO 



n=0 



(21) 



Note that K-matrix is a generalization of an idea first pro- 
posed by Champagne and Kirtman.^^ If the sums rules 
are exactly obeyed, the K-matrix vanishes. Therefore, 
kappa is calculated for a given system by substituting 
the known values of the transition moments and energies 
into Equation 1211 Deviations of Kmp from zero is a sign 
that either the act of truncation in inaccurate; or that 
the matrix elements of the dipole operator and energies 
are incorrect for that particular system. 

In order to test the pathologies associated with trunca- 
tion, we must start with a potential energy function for 
which the Schrodinger Equation can be solved analyti- 
cally. This approach ensures that the energies and dipole 
matrix elements are physically sound and that patholo- 
gies or inaccuracies inherent in approximation techniques 
are avoided. The clipped harmonic oscillator (CHO) 
is the ideal system (where the potential is harmonic for 
X > and infinite for x < 0) since it is the simplest 
case of an asymmetric potential that yields a hyperpo- 
larizability that is near the fundamental limit. 10] The 
matrix elements of the position operator of the clipped 
harmonic oscillator (CHO) are given by, 



_ ^MAX 
•^mn — '''10 yran-! 



(22) 



where the dimensionless matrix gmn is defined by 



Smn — I — ( 1) 



((m+n)/2) 



7,1 



(m-n)2-iy \^/m\n\, 
(23) 



where to and n are restricted to the odd integers. The 
energy for state n is given by 



En = hWQ 



Figure El shows a surface plot of Kmp 



(24) 



and the inset 

shows the matrix elements. Recall that the K-matrix is a 
measure of the fractional deviation of the particular trun- 
cated sum- rule equation Xnm- As such, when truncating 
to three levels, the sum rule A'qq with the clipped har- 
monic oscillator wavefunctions gives a deviation of less 
than 7% from the infinite-level case. The sum rule equa- 



tions x'"y^ and X22 , which we dropped from our analysis, 
deviate by 82% and almost 270%, respectively. 

Since the clipped harmonic oscillator's hyperpolariz- 
abihty is near the fundamental limit, we would expect 
that the three-level sum rules that are kept in the anal- 
ysis to deviate substantially less than the equations that 
are dropped. Figure |3| confirms that the domain of the 
sum rules used in our maximization calculation is rela- 
tively flat. 

To show how quickly the sum rules converge, and how 
well the exact clipped harmonic oscillator obeys them, we 
consider a 32-level model. Figure O shows a contour plot 



.(3) 



of K 



(32) 
nip 



and 



J32) 
^rap 



First, we note that the deviations 



alternate sign when m + p changes from even to odd. 

(3) 

(This same pattern is also observed for K,np in Figure 
121) The jagged line in Figure represents the point at 



which K' 



(32) 
mp 



0.025. The areas that are shaded black 



represented matrix elements of n that are larger than 0.1. 
As such, most of the sum rules are accurate. Indeed, even 



.(10) 



converges quickly and the deviations are small for 
all sum rule equations for m + p < 10. 

Figure|3|D shows the absolute magnitude of 
more clearly shows deviations of the K-matrix. In the 



^,(32) 



which 
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that the K-matrix is far larger than one would expect if 
it depended linearly on the random error. In the white 
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FIG. 3: Contour plot of (a) Kmp and (b) 



white region, > 0.1. Above the diagonal, from 

(m,p) = (1,30) to (m,p) = (30,1), n'^p' > 0.03. Re- 
call that in the case of the three-level model, we ignored 

(3) . . (3) 

the two equations: for contradicting Xqq and for 

f31 

being unphysical; and Equation since it yielded un- 
physical results. The K-matrix may be a useful tool in 
studying which sum rule equations that need to be dis- 
carded in higher-level models. 

It is interesting to speculate whether the abrupt sign 
changes in the K-matrix may result in a cancellation of 
terms in the SOS expression for the nonlinear-optical sus- 
ceptibility; and, if small errors in the matrix elements 
may lead to large uncertainties in the nonlinear-optical 
response. 

To test this hypothesis, the matrix elements and ener- 
gies of the clipped harmonic oscillator were recalculated 
with random fluctuations of ±5 % added to each using 
the equations. 



and 



En 



E„ 



(1 + r), 



(25) 



(26) 



where r is a random number equally distributed in the 
interval —0.05 < r < -1-0.05. Figure^lshows a plot of the 
magnitude of the K-matrix for these parameters. Note 




FIG. 4: Contour plot of jftmp'! where dipole moments and 
energies randomly are changed by ±5%. 



region, Kmp^ > 0.083, while in many regions, k' 



(32) 
nip 



> 1. 

So, it appears that a small error in the matrix elements 
and energies can lead to a large violation of the sum rules. 



IV. CONCLUSION 

In conclusion, we have shown that the sum rules 
can be used to calculate the fundamental limits of the 
nonlinear-optical susceptibility provided that the system 
approaches a three-level model as (3 and 7 approach the 
fundamental limit. This process of truncating the sum 
rules and demanding that the fundamental limit occurs 
in the three-level limit is motivated by exact results that 
are calculated for the polarizability. We introduce three 
guiding principles that aid us in determining which sum 
rule equations should be dropped, and show that the fun- 
damental limit is unique. Furthermore, we introduce a 
method for testing the validity of the transition moments 
and energy levels and find that even small errors in them 
can lead to large violations of the sum rules. As such, this 
method is ideal for testing the results of semi-empirical 
calculations. 
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